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Abstract 



Motivated by numerous X-ray scattering studies of lamellar phases with 
membrane proteins, amphiphilic peptides, polymers, or other inclusions, we 
have determined the modifications of the classical Caille law for a smectic 
phase as a function of the nature and concentration of inclusions added to it. 
Besides a fundamental interest on the behavior of fluctuating systems with 
inclusions, a precise characterization of the action of a given protein on a lipid 
membrane (anchoring, swelling, stiffening ...) is of direct biological interest 
and could be probed by way of X-ray measurements. As a first step we con- 
sider three different couplings involving local pinching (or swelling), stiffening 
or tilt of the membrane. In the first two cases we predict that independent 
inclusions induce a simple renormalization of the bending and compression 
modulii of the smectic phase. The X-ray experiments may also be used to 
probe correlations between inclusions. Finally we show that asymmetric cou- 
pling (such as a local tilt of the membrane) results in a modification of the 
usual Caille law. 
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Lyotropic smectic phases (L a phases) are liquid crystalline systems of parallel, regularly 
spaced lipid bilayers Those phases form a quasi-crystalline structure in one (the z) 

direction, while retaining their fluid properties in the two others directions (x — y plane). 
They have been studied extensively in recent years, as they form one of the most convenient 
lyotropic structure for the experimental study of fluctuating interfaces. They have been for 
instance particularly useful in demonstrating the importance of the so-called entropic (Hel- 
frich) repulsion between fluctuating membrane [Q. Furthermore, the degree of dilution (the 
layer spacing) and other controlled parameters of the lyotropic L a phase can be experimen- 
tally varied over a wide range, which allows for a precise study of the scattering intensity 
of the phase. This scattering intensity is expected to show several peaks corresponding to 
a stack of regularly spaced membrane. In these layered systems however, the thermal fluc- 
tuations destroy the long range order at finite temperature, leading to weak - power law - 
singularities of the peaks. 

Inclusions in membrane have been actively studied for many years. Most of the earlier 
work was done by the biological community, for which it was a necessary step toward the 
comprehension of the complexity of biological (cell) membrane In the last twenty years 
or so, a more physical description of membrane-inclusion complexes, in term of curvature 
energy and disruption of the molecular order of the membrane, have led to interesting results 
on membrane-induced interactions between inclusions ||. While this "physical" picture of 
the membrane-inclusion interactions leaves aside many fundamental processes driven by 
specific interactions, it allows for a description of universal features of the complexes in 
terms of elastic couplings between membrane and inclusion. This description is applicable 
to a certain extent to membrane with adsorbed or grafted polymers or colloids H0], as well 
as to transmembrane inclusions resembling certain membrane proteins ||. 

In recent years, numerous experimental studies have investigated the role of added mem- 
brane proteins in lyotropic smectic phase. This lipid phase seems to be of particular interest 
for the study "in vitro" of the bilayer-protein association P-TT2]], and may have some prac- 
tical application for gene therapy ||13|| . From a physicist point of view, it is interesting to 
understand the effect of inclusions on the thermal fluctuations of flexible, fluid membranes. 
From a more "biological" point of view, a precise characterization of the physical action of 
a given protein onto a lipid membrane (anchoring, swelling, stiffening ...) is of direct biolog- 
ical interest, and could be reach via X-ray measurements. This work have some relevance 
in material science as well, since for practical purposes, a lamellar phase is very rarely com- 
posed of two components (solvent -l-surfactants) only. The adjunction of a third or fourth 
component, such as cosurfactant, polymers, proteins or colloids is of common practice to 
tune the properties of the phase [|TJ,|rj]. The object of the present work is to determined the 
modifications of the smectic X-ray scattering due to the presence of inclusions in the lamellar 
phase. Ultimately, one could expect to map out different membrane-inclusion couplings by 
their effect on the X-ray structure factor. We consider three model inclusions, leading to 
different modifications of the membrane elastic properties, and we show that the qualitative 
analysis of the diffraction spectrum is enough to discriminate between these inclusions. 

Several recent theoretical works deal with lamellar systems with flexible polymers |15 



or other soluble inclusions (so-called doped solvent lamellar phases). These works are com- 
plementary, as we are interested in insoluble inclusions such as the membrane proteins used 
in H, while [[TJJ studied effects such as polymer depletion - induced interactions between 
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lamellae, resulting from an equilibrium between the trapped inclusions and free inclusions 
in pure solvent. 



^> The standard description of Caille [Hi for pure L a phases has been confirmed experi- 



mentally in numerous systems. It is based on the classical smectic hamiltonian fl7|| : 



(1) 



where u(r) is the (continuous) normal displacement of the layer at the point r, z is the 
coordinate normal to the layers, and Vj_ is the gradient along the layers. The elasticity of 
the smectic phase is characterized by the two coefficients K and B, which are related to the 
energy cost of bending and compressing the sample. Typically, K = n/d, where k is the 
bending modulus of a single membrane and d is the average layer spacing, and B is function 
of the intermembrane interactions Those modulii define the characteristic smectic length 
||17|| . The scattering intensity /(q) is related to the Fourier transform of the 



A = y/K/B 

density correlation function G n (r), itself related to the exponential of the layer displacement 
correlation function g(r) = go/2((w(r) — w(0)) 2 ), where (...) represents a thermal average of 
the fluctuations, and go = 2ir/d (see @, sec 6.3.2). As usual, the Fourier decomposition of 
the layer displacement is used: w(r) = / c? 3 g/(27r) 3 M g e* qr . Using Eq.(ffl): 
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where V is the volume of the sample. After integration over q z , an asymptotic solution of 
g(r) can be found [18|: 
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where Jo(x) is the Bessel function of the first kind, 7 is the Euler constant, and E\(x) is the 
exponential integral (/ is a molecular size). From Eq.(||), one can show that the scattering 
intensity I n (q) near the n th peak, defined by 



I n (q) = G n (q + nq ) + G„(q - nq ) with G„(r) oc e 



-n 2 g(r) 



shows power law singularities of the form p~6|JT9[] 
G„(q- nq ) - 
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(5) 



The above law gives a satisfactory description of the X-ray structure factor in the vicin- 
ity of the diffraction peaks, and has been observed experimentally (see |I| and references 
therein). A more accurate calculation of /(q) over a broad range of wave vector requires 
a discrete description of the lamellar phase [p0|1 , and the account of the fluctuations of the 
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lipid concentration and is not tractable analytically. As one can see, this peculiar power 
law divergence allows for an experimental determination of the product KB, through the 
parameter 77. Low values of the parameter 77, such as those observed for electrostatically 
stabilized systems (~ 0.2) correspond to a well defined smectic organization with little fluc- 
tuations. Larger values of order unity are observed for weakly interacting neutral or screened 
systems stabilized by the Helfrich interaction fLf. 



<0> Several couplings between the inclusions and the membrane are considered below. 
They are referred to as "pinch", "stiff", and "tilt", and are schematically depicted in Fig.l 
and Fig. 2. The "pinch" inclusion can be thought of as exerting a force which tends to 
pinch (Ppi nc h > 0, see Eq.(|6]) below) or swell (/3 p i nc h < 0) neighboring membranes. Since we 
are using a continuous description for the lamellar phase, the inclusion may reside within 
a membrane or between two membranes. Inclusions which would typically induce such 
deformations are amphiphilic proteins [|§, with hydrophobic parts which penetrate in the 
bilayers and hydrophilic parts lying in the water, or simply particles larger than the layer 
spacing which sterically swell the membrane. A pinched lamellar structure has also been 
identified in the case of anionic polymers added to a lamellar phase formed by a mixture of 
cationic and neutral lipids |22| . The hamiltonian corresponding to this perturbation couples 



the density of inclusion p(r) with the variation of the layer spacing ||23|| : 

AH pmch = Jd^ pmchP (r)d z u (6) 

The "stiff" inclusion corresponds to a transmembrane protein which locally stiffens (or 
softens) the membrane and modifies its fluctuations. A stiffening is expected for many 
integral membrane proteins, or for rigid inclusions laying on the lipid bilayer. On the other 
hand, pores in the membrane are expected to induce a softening of the membrane (among 
other, less trivial phenomena). The inclusion density is coupled with a variation (SK) of the 
bending modulus p4| , |25|| : 

AH sttff = ^J d 3 r6Kvp(r) (v 2 ± u) 2 (7) 

where v is the effective volume of an inclusion in the membrane. 

The "tilt" inclusion corresponds to an anisotropic inclusion (such as a grafted polymer) 
which would locally induce a curvature to a membrane without changing the layer spacing. 
It has been shown that such a deformation can be described by a modification of the "pinch" 
hamiltonian Eq.([|) p6| ; we will get back to it later. It is clear that any real inclusion will, 
to a certain extent, combine the three effects described above. 

A single inclusion of the "pinch" or "tilt" type induces a non-zero average deformation 
to the membrane. However, this mean deformation evens out during the course of an 
experiment, and is not detectable on the scattering spectrum (except if it induces a phase 



separation, such as expulsion of the solvent [11,23]). The effect of the inclusions on the 



spectrum is visible via their effect on the membrane fluctuations. In what follows, we 
do not investigate a possible coupling between the membrane fluctuations and the inclusion 
correlation function (assumption of a quenched distribution of inclusions). It has been shown 
for a single membrane system that "quenched" and "annealed" distributions of inclusions 
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give similar corrections to the fluctuation spectrum of the membrane [25]. The scattering 



experiments show an average of the structure factor over the positions of the inclusions, it 
is thus necessary to calculate the membrane correlation function after a proper average over 
the possible positions of the inclusions. 



<0> For the case of "pinch" inclusions, the hamiltonian Eqs. 
discrete form: 



can be rewritten in the 
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The structure factor is then easily obtained: 



<ki 2 > t , uyy 



V 



Q q 



+ 



Qo 



k B T 2 q 2 z {\ Pq \ 2 ) f 

Bql + Kqi +P ^ c \ B ql + Kq\) 2 



V 



(8) 



(9) 



where {■■■) p designates an average over the position of the inclusions. The first part of the 
left hand side corresponds to the correlation function in a pure L a phase ^°)(r) (Eq.©), 
and the second part describes the corrections due to the inclusions Ag. This correction 
reflects correlations between inclusions (|p g | 2 ) p - 

If the inclusions are independent from each other, (\p q \ 2 )/V is the average concentration 
of inclusions p. Similarly to Eq.(§), the correction to the spatial structure factor can be 
written: 
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This correction leads to the same power law divergence of the scattering intensity as the 
pure L a (Eq(H)), with the effective modulii: 
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One notice right away that the "pinch" inclusions tend to soften the lamellar phase. We 
thus expect broader peaks in the scattering intensity (see Fig.l). This finding is qualitatively 
consistent with experiments on comparable systems Note that results reported for 

the doped solvent situation tend to show qualitatively similar results for physically different 
reasons. In that case, the inclusions-mediated attraction between lamellae tends to reduce 
the compression modulus B [ |T5|| . 

Possible correlations between inclusions may be induced either by direct interactions 
between them, or by membrane mediated interactions. In the Latter case, it has been 
shown that the inclusions interact via the potential (in Fourier space) |23j 
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In this case, a "Debye-Huckel like" j| expansion of the correlation between inclusions (valid 
if pU q /{k B T) <C 1), shows that: 



(IPgl 2 ) P Qr (\P,\ 2 ) = PQq 



V 1 + pU q / (k B T) V Q' q 

with Q q ,=B'ql + Kqi and B' = B(l-^^j (13) 

From Eq.(|9p, the corrections to the structure factor is given by: 

A(K| 2 > = %^ - r)pinc h = v(l+ , ~ 3^ n*\ ( 14 ) 
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The effective modulii are now given by: 

KB 

( KB ) P inch = -/ - : — — w ( 15 ) 

\ 2k B TB ~r 16 I 2kgTB J J 

Note that in the limit of validity of the expansion, the correction due to correlation is quite 
small. 



<0 The investigation of the effect of the "stiff' inclusions is more delicate, because of 
the non-quadratic nature of the hamiltonian Eq.(^). This leads to a coupling between 
different Fourier modes of the membrane displacement, in which case Eq.(§) for the real 
space correlation function has to be replaced by: 

9(r) = \* I Wfw?^ (e ' ,r " ( e "' r " 1} (16) 

The discrete version of the hamiltonian reads: 



H. 
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= i £w (vSr«W + ^ff^iV^o) UqU « (17) 

from which the mean square displacement can be calculated: 
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(18) 



In the limit where the correction is only a small perturbation, an expansion of the previous 
equation [G° + AG}^, = [G%^ - [AG] q>q , x [G% 2 q , leads to: 

(u q u*,) _ k B T 5Kvk B T q±q'±Pq-g> nQ x 



V Bq 2 z + Kqi q ' q V (Bq 2 z + Kq\_) (Bq' z 2 + Kq' ± 4 ) 

For an homogeneous distribution of particle of average concentration p, the Fourier modes 
are decoupled: p q = p5 q . It is then easy to show that again, the scattering intensity show 
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peaks with power law divergence, from which the product of the effective compression and 
bending modulii can be extracted. The renormalization introduced by the inclusions is given 
by: 

(KB) «'» = z He? (20) 

Note that this results compares to the finding of Netz and Pincus |]25|| , who have shown 
that for a single membrane with inclusions, the effective bending rigidity of the membrane 
is given by K e ff ~ k + 8k(/) at the lowest order in the volume fraction <fi of inclusions. 

In contrary to the case of "pinch" inclusions, a stiffening of the L a phase is observed 
for inclusions increasing locally the bending rigidity of the membrane. This should lead to 
a sharpening of the peaks in the scattering intensity. Hence, a quick look at the intensity 
scattered by the sample should give us valuable information on the inclusion-membrane 
coupling. A stiffening of a L a phase has been recently observed upon adjunction of polysoap 



molecules [pq| . This is a more complex case, since the polysoaps are fairly long molecules 
and can be expected to modify (increase) the interactions between membranes in addition 
to their effect on an isolated bilayers. The reduction of the effective r] has however been 
observed below the overlap concentration of polysoaps, where the molecules are expected to 
be spread flat onto the membrane. 



^> The last coupling investigated in this work is the case of anisotropic inclusions which 
induce a local tilt (spontaneous curvature) on a membrane. It has been shown in a previous 
work p6[ that such a tilt can be conveniently mimic with the aid of an analogy between 
the "pinch/swell" inclusions and +/- electrical charges in electrostatics (see Fig. 2). The 
deformation of the smectic induced by a "dipole" p z = (3 t utd (where j3 t ut, the tilt intensity, 
has the dimension of an energy) normal to the membrane has been calculated [26|: 



and correspond to imposing a local curvature C = /3 t ut/ ij^d) at the position of the inclusion 
(r_i_ = 0, z = d/2) without changing the layer spacing (k = Kd is the rigidity of an isolated 
membrane). It can be shown easily in the Fourier space that such a deformation can be 
induced by a coupling between the second order z-derivative of the layer deformation and 
the inclusion density: 

AH tilt = J d 3 r PzP (v)d 2 z u (22) 

which, for independent inclusions, leads to the membrane structure factor:: 

(Kl 2 ) k B T qt 

V ~ B q l + K q i +Pz {Bq l + K q if P {I6) 

and the correlation function: 
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After integration over q z (the (_L) indexes are forgotten): 



where J = 3 / c% 3 (l - J (gr)e^ W ) + Xz J dqq b J Q {qr)e~ Xzq2 (25) 



Both integrals can be computed after a series expansion of the Bessel function. The correc- 
tion to the correlation function takes the form: 

A ^ (r) = ( j^f ( x + u ~ u2 ) + constant ) with u = ^ (26) 

The X-ray spectrum is obtained by a combination of the above expression with Eqs.(|3|j4]): 
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Unlike the case of symmetrical inclusions considered earlier, the asymptotic scaling laws 
for G q in the limits q z = or qj_ = cannot be easily extracted from the above expression, 
since the exponential term in the integral cannot be written as a function of r 2 ± /{Xz) alone. 
However, up to the first order in the dipole strength, Uat/z < 1 for z > d/2 and we can 
expand the exponential in the integral Eq . (]27|) . One can see clearly that unlike symmetrical 
inclusions such as "pinch" and "stiff", the effect of a "tilt" inclusion cannot be matched 
by a renormalization of the smectic modulii, as they modify of the power law divergence of 
the scattering intensity near the Bragg positions. In fact, the inclusions add a term to the 
intensity, which for q± = can be written: 

W)oc-g2 4/^/^e^/^e-(l + a;-a; 2 ) e -^ ; Z es q z z (28) 

and in the other limit q z = 0: 

W^)«-^4/^/^Jo(^)/^e--(l + a;-a; 2 )e-^H ; R = q ± r ± (29) 

The integrals being prefactors only, the scattering intensity should take the form: 

T nr / a ^7 2+v ~ a 2<f z for Q± = (o n \ 

J-tilt \ -4+2r? 2n r n [<JU 

{ aig ± - a 3 q ± ' for q z = 

where the 04 are constants. The first term is the Caille law for a pure smectic phase and 
the second term is the signature of the asymmetric inclusions. Note that the effect of the 
tilt inclusion as described by Eq.(^) is most probably very small and might be hardly 
detectable in experiments. Note also that as a "pinch" inclusion lead to a renormalization 
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of the compression modulus B (see Eq.([T^)), a "tilt" inclusion can be seen as renormalizing 
the elastic modulus of an higher order term of the smectic hamiltonian, of the form K'(d 2 u) 2 , 
which has been neglected in the reference hamiltonian Eq.fll]) P9 |. This reinforces the fact 



that their effect on the structure factor must be small. One nevertheless might hope to 
detect such an effect because of its dependence on the concentration of inclusion. 

<0> In conclusion, we have derived the effect of inclusions on the shape of the scattering 
curve of a smectic L a phase. We have considered the three inclusions depicted in Figs.l and 
2, and have shown that the two symmetrical inclusions ("pinch" and "stiff") renormalize 
the elastic modulii of the phase without modifying the power law divergence near the peaks 
(Caille law). In agreement with experiments, we have obtained that the "pinch" inclusions 



soften the lamellar phase flllj], while the "stiff" inclusions makes it more rigid fL2| . the "tilt" 
inclusions are predicted to modify the shape of the scattering curves, but their effect is small 
and might hardly be detectable. 

As a last remark, we note that in view of our results, the study of the X-ray scattering 
spectrum of a lamellar phase with inclusions should be a powerful tool to investigate some 
changes in conformation of membrane proteins. One can argue that the bridge (pinch) 
between neighboring membranes created by some amphiphilic peptides and reported in pTj 



may be pulled out of one of the two membranes as the intermembrane spacing is increased. 



If the peptide is flexible enough, it is then expected to lay flat on one membrane fll2 |. This 
conformation change represents a transition from a "pinch" to a "stiff" inclusion, and should 
be directly detectable from the shape of the scattering peak. Moreover, some amphiphilic 
helical peptides in contact with a lipid bilayer are known to go from an adsorbed state at 
low peptide concentration, to an inserted state (a certain number of peptide aggregating to 
form an hydrophilic pore) at higher concentration |30| . This transition from one state to the 
other could be interpreted in terms of "stiff" and "soft" (in the case of a pore) inclusions; 
One could investigate X-ray scattering from such a solution with this transition in mind. 
Nevertheless, one should always bear in mind that the different states of the amphiphilic 
peptides are most probably a complex mixed between the types of inclusion proposed in our 
work, and probably others. 
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Fig.l The two symmetrical inclusions considered in the paper. A "pinch" inclusion modifies 
the layer spacing and broadens the Bragg peak - A "stiff" inclusion increase the local rigidity 

and makes the Bragg peak narrower 
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(a) 



(b) 



Fig. 2 a) The "tilt" inclusion induces a spontaneous curvature 
b) How to make a "tilt" from the combination of "pinch" and "swell 
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